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INTRODUCTION 


Recent  earthquakes,  particularly  those  in  Alaska,  Japan,  and  Chile, 
have  emphasized  the  high  damage  threat  the  soil  liquefaction  phenomenon 
poses  to  waterfront  structures.  These  experiences  have  shown  that  both 
the  nature  of  waterfront  facilities,  such  as  earth-retaining  structures, 
and  the  depositional  environment  of  the  coastal  marine  soil  contribute 
to  major  liquefaction  damage.  A  recent  study  conducted  by  the  Office  of 
Naval  Research  (Ref  1)  recognized  a  major  liquefaction  hazard  existing 
at  West  Coast  Naval  Stations.  A  more  recent  investigation  at  the  Naval 
Air  Station  (NAS)  North  Island,  Calif.  (Ref  2),  concluded  that  liquefac¬ 
tion  under  design  earthquake  levels  could  result  in  destruction  of  such 
critical  structures  as  aircraft  carrier  berths,  aviation  fuel  tank 
farms,  and  underground  utility  service  lines.  Unfortunately,  almost  all 
previous  studies  of  the  liquefaction  problem  have  been  concerned  with 
either  conventional  building  foundations  or  with  analyses  of  dams,  and 
extensive  information  is  not  available  on  the  effect  of  liquefaction  on 
these  types  of  specialized  structures. 

Basically  two  approaches  are  used  for  evaluating  the  liquefaction 
potential  of  a  deposit  of  saturated  sand  subjected  to  earthquake  shaking 
(Ref  3). 

1 .  Field  Data  of  Performance  of  Sand  Deposits  in  Previous  Earth¬ 
quakes  .  Post  earthquake  surveys  of  areas  where  liquefaction  has  or  has 
not  occurred  have  been  used  to  prepare  charts,  based  primarily  on  the 
Standard  Penetration  Resistance  of  the  deposit,  for  differentiating 
between  liquefiable  and  nonliquef iable  conditions.  Empirical  comparisons 
and  evaluations  of  this  type  do  not  take  account  of  such  significant 
factors  as  the  duration  of  shaking  or  the  extent  of  drainage  and  depend 
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upon  the  reliability  of  field  observations  and  such  field  tests  as 
penetration  resistance  (often  after  the  fact).  Thus,  many  engineers 
feel  that  such  correlations  provide  only  preliminary  evaluations  of 
liquefaction  potential.  These,  they  feel,  will  often  need  to  be  supple 
mented  by  detailed  studies  based  on  ground  stress  analyses  and  soil 
testing  programs. 


2.  Evaluation  of  Stress  Conditions  in  the  Field  and  Determination 


of  the  Stress  Conditions  Causing  Liquefaction  of  Soils.  Analytical 
procedures  for  evaluating  the  liquefaction  potential  of  soil  deposits 
involve  two  independent  determinations:  (a)  an  evaluation  of  the  cyclic 
stresses  induced  at  different  levels  in  the  deposit  by  the  earthquake 
shaking  and  (b)  investigation  to  determine  the  cyclic  stresses  that,  for 
given  confining  pressures  representative  of  specific  depths  in  the 
deposit,  will  cause  the  soil  to  liquefy  or  undergo  various  degrees  of 
cyclic  strain.  The  evaluation  of  liquefaction  potential  is  then  based 
on  a  comparison  of  the  cyclic  stresses  induced  in  the  field  with  the 
stresses  required  to  cause  liquefaction  -  or  an  unacceptable  limit  of 
cyclic  strain  -  in  representative  samples  in  the  laboratory.  This 
approach  raises  questions  regarding  the  validity  of  the  laboratory 
results  with  respect  to  field  conditions  and  is  essentially  limited  to 
level  ground  situations  away  from  load  discontinuities. 

The  Navy  has  many  facilities  on  reclaimed  or  marginal  land  in 
waterfront  locations.  The  types  of  structures  the  Navy  builds  on  such 
lands  (e.g.,  quay  walls,  drydocks ,  storage  tanks,  etc.)  are  not  widely 
studied  in  the  civilian  sector.  Current  analytical  techniques  are  not 
adequate  regarding  analysis  of  conditions  other  than  those  involving 
level  ground  away  from  structures.  The  uncertainty  associated  with  soil 
strength  ami  earthquake  occurrence  is  not  routinely  used  to  estimate 
liquefaction  potential.  It  is  apparent  that  the  precision  with  which 
liquefaction  investigations  can  be  carried  out,  particularly  with  regard 
to  ocean  front  facilities,  leaves  much  to  be  desired.  Better  evaluation 
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t  echniques  must  In'  developed  for  determi  it  1  it#  !  i  que t  act  ion  hazard  m  the 
mailin'  depos i t i ona l  environment  i-xist  ui|!  at  many  Naval  est  al>l  i  shnii'ut  s  . 

Out’  ot  tin-  majoi  short  com  mgs  m  av.iil.ih  It*  liquefaction  evaluation 
technology  is  related  to  the  lack  ot  validated  or,  in  some  cases,  even 
credible  procedures.  It  is  apparent  that  realistic  large-scale  tests 
for  investigating  the  liquefaction  phenomenon  are  desirable.  However, 
because  ot  the  greater  number  ot  complex  variables  which  control  lique¬ 
faction  potential,  the'  results  to  be  derived  from  an  individual  large- 
scale  lest  at  this  time  might  not  warrant  the  very  large  expenditures 
that  the  Navy  would  incur.  I'tus  study,  therefore,  is  currently  examining 
the  overall  field  of  liquefaction  prediction  and  evaluation  and  concen¬ 
trating  upon  a  broad  program  on  the  feasibility  ot  adapting  available 
and  anticipated  liquefaction  analysis  developments  to  Navy  requirements. 

1'his  study  concentrates  on  two  main  aspects  ot  liquefaction  research: 
ill  better  site  evaluation  technology  and  12)  realistic  liquefaction 
analytical  techniques.  Earthquake  prediction  and  ground  motion  studies 
are  outside  the  scope  of  this  project. 

NAS  North  Island,  which  had  recently  been  designated  as  having  a 
potential  liquefaction  hazard  (Net  2),  was  selected  for  furthei  study. 
Several  problem  soils  at  this  facility  were  sampled  using  fixed  piston 
sampling  techniques;  the  most  critical  soil  was  subjected  to  cyclic 
triaxial  testing.  flic  results  ot  this  testing  procedure  are  discussed 
m  the  next  section  ot  this  document.  fhe  North  Island  site  will  be 
used  toi  further  evaluation  ot  soil  reconnaissance  techniques  applied  to 
identifying  liquefaction  potential.  In  order  to  carry  out  this  program 
more  effectively,  CE1.  is  now  able  to  conduct  cyclic  triaxial  tests, 
fins  will  permit  in-depth  evaluations  ot  the  liquefaction  potential  ot 
characteristic  waterfront  soil  deposits  using  the  current  technology. 

Analytical  procedures  tot  predicting  liquefaction  response  have 
been  improved  by  use  of  effective-stress  computer  code  analysis  Compute i 
codes  are  discussed  further  in  a  later  section  ot  this  document.  A 
state-of-the-art  guide  (.Net  A)  dealing  with  structures  located  in 
ea i thquake-prone  areas  has  been  published  tor  use  by  the  Naval  Shore 
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F.s  t  .»t>  1 1  s  hiniMt  t  .  1  In  s  guide  Givi's  general  background  data  on  the  lique¬ 

faction  problem  and  provides  .1  summary  of  current  technology.  Reference 
■4  was  excerpted  largely  from  a  more  comprehensive  series  of  reports 
completed  tot  the  Federal  Highway  Administration  (Ret  5-7). 

File  objective  ot  the  research  program  discussed  herein  is  to  develop 
procedures  tor  reliably  and  quant  1 1 at i vely  evaluating  the  threat  levels 
posed  by  earthquake-generated  soil  l 1 que t act i on  at  waterfront  facilities. 
These  procedures  and  any  associated  instrumentation  development  will  be 
suitable  for  1 ncorpora t i on  into  a  Navy  manual  tot  designing  new  facilities 
and  taking  t  hreat -1111 1  i  gat  1  ng  measures  tor  existing  ones. 


NORTH  l SLANT)  SOILS 

Outing  the  investigation  ot  the  earthquake-induced  liquefaction 
threat  at  NAS  North  Island  tRot  21,  several  uncommon  soils  were  encoun¬ 
tered.  Port  10ns  of  t  In'  Air  Station  originally  below  sea  level  had  been 
reclaimed,  using  hydraulic  tills  -  primarily  siltv  sands  -  dredged  from 
the  bottom  ot  San  Diego  Bay.  These  reclaimed  areas  had  been  subjected 
to  several  series  ot  tilling,  most  between  1'HU  and  ldS2.  As  a  result 
ot  invest igat ions  made  on  North  Island  in  connection  with  foundations  ot 
existing  structures,  a  recently  tilled  region  was  selected  foi  further 
study  (Ret  2).  Available  soil  data  were  supplemented  by  additional  soil 
investigations,  including  dynamic  split  spoon  penetration  tests  and 
statu  cone  penetrometer  soundings.  The  split  spoon  penet 1  at  1 ons  were 
conducted  as  specitied  111  ASTM  1)  lf>So  except  that  .1  2-1;  2- inch  (t'a  1  i  t orni .1- 
type)  split  spoon  was  used.  The  cone  penetration  soundings  conformed  to 
ASTM  D  .<•♦*«  1-7.ST.  Based  upon  this  111form.it  ion  the  generalized  soil 
protile  shown  111  Figure  1  was  developed  Tins  soil  profile  encompasses 
an  old  tilled  ba\  channel,  formerly  known  as  Spanish  Bight,  which  was 
hydraulically  tilled  to  its  present  elevation  111  l')4!S.  The  old  bnv 
bottom  immediately  prior  to  the  tirst  tilling  is  represented  by  layer  .S . 
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I'he  results  of  split  spoon  penetration  tests  am)  friction  cone 
soundings,  denoted  as  penetration  holes  PI  through  PS  along  the  soil 
profile  of  Figure  1,  are  shown  in  Figure  2.  It  should  be  noted  that  the 
split  spoon  penetration  values  near  the  boundary  between  layers  4  and  S 
tall  as  low  as  one  or  two  blows  per  foot.  Figure  3  (Ref  3)  shows  plots 
of  standard  penetration  test  results  (number  of  standard  hammer  blows 
required  to  drive  a  standard  sampler  12  inches)  versus  in  situ  soil 
relative  density.  For  a  granular  soil  at  confining  pressures  of  5  to  10 
psi,  Figure  3  suggests  a  relative  density,*  D  ,  for  soil  stratum  5  of 
less  than  50%. 

In  Figure  4,  Schmertmann  (Ref  9)  presents  relationships  between 
relative  density  and  static  cone  penetration  resistance.  With  the 
relatively  low  friction  cone  penetration  readings  of  5  to  10  kg/sq  cm 
shown  in  Figure  2,  values  of  relative  density  on  the  order  of  only  20% 
or  30%  are  indicated  for  stratum  5. 

As  suggested  in  Reference  2,  natural  soils  with  relative  densities 
of  less  than  40%  are  expected  to  be  quite  rare  in  nature.  However, 
based  upon  the  available  information  it  was  necessary  to  assume,  for 
analysis  purposes,  that  the  soil  in  the  region  of  layer  5  could  have  a 
relative  density  no  greater  than  35%. 

The  resistance  of  granular  materials  to  cyclically  induced  liquefac¬ 
tion  is  a  function  of  the  in  situ  relative  density.  Relative  densities 
less  than  about  45%  are  generally  considered  quite  critical  and  can 
permit  disastrous  liquefaction  failures  at  fairly  low  ground  acceleration 
levels.  For  example.  Figure  5  shows  a  record  of  previous  liquefaction 
occurrences  as  a  function  of  horizontal  shear  stress  over  confining 
stress  ratio  and  relative  density  (Ref  10).  This  figure  suggests  that 
for  relative  densities  on  the  order  of  40%,  stress  ratios  of  about  0.1 
would  probably  result  in  liquefaction.  A  stress  ratio  of  0.1  corresponds 
to  horizontal  ground  accelerations  of  about  0.08  (Ref  11)  times  the 
acceleration  of  gravity  (or  modified  Mecca lli  intensities  of  VI  or  VII). 

'■'Relative  density  is  a  measure  of  the  present  density  of  a  soil  with 
respect  to  its  maximum  and  minimum  densities. 
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Since  most  regions  in  California  can  be  expected  to  experience  ground 
motions  of  this  severity  relatively  often  (in  the  geologic  time  scale), 
it  was  considered  critical  that  the  liquefaction  potential  of  this  soil 
be  explored  further.  This  situation  is  further  amplified  when  one 
realizes  that  this  type  of  coastal  soil  commonly  referred  to  as  "bay 
deposits"  seems  to  be  fairly  widely  distributed  along  the  coast  of 
Southern  California  and  appears  at  other  Naval  facilities  such  as  Long 
Beach  Naval  Shipyard. 

Therefore,  a  program  of  undisturbed  sampling  and  laboratory  cyclic 
triaxial  testing  was  carried  out.  The  sampling  program  consisted  of 
obtaining  soil  samples  along  the  soil  profile  shown  in  Figure  1  using  an 
Osterberg  piston  sampler  (Ref  12).  This  sampler  utilizes  a  thin-walled 
sampling  tube  3  inches  in  outside  diameter,  2.88  inches  in  inside  diameter, 
and  3b  inches  in  length.  A  complete  description  of  this  sampling  and 
testing  program  is  presented  in  Reference  13.  The  undisturbed  samples 
were  handled  and  transported  to  the  laboratory  using  approved  techniques. 
From  these  samples,  characteristic  specimens  were  selected  for  cyclic 
triaxial  testing. 

The  soil  specimens  were  subjected  to  the  same  level  of  confinement 
as  the  estimated  initial  in  situ  vertical  effective  stress,  and  then 
saturated  to  "B-values"  greater  than  0.97  prior  to  testing.  Figures  6 
and  7  present  the  results  of  cyclic  triaxial  testing  on  eight  soil 
specimens  extracted  from  soil  layer  5.  Two  criteria  for  soil  failure 
were  used.  Figure  6  presents  maximum  shear  stress  to  confining  stress 
ratio  versus  the  number  of  load  cycles  that  result  in  an  axial  or  vertical 
strain  of  5%.  Figure  7  presents  the  same  type  of  information  but  is 
based  upon  the  time  at  which  measured  pore  water  pressure  in  the  specimen 
equals  the  triaxial  confining  or  chamber  pressure  (rather  than  the  5% 
strain  level).  Both  criteria  lead  to  the  same  result  in  this  case. 

In  Figure  8  are  shown  the  North  Island  cyclic  triaxial  data  for 
soil  layer  5  plotted  on  top  of  shake  table  data  presented  in  Reference 
10.  The  triaxial  data  have  been  plotted  in  terms  of  the  same  stress 
ratio  used  for  the  shake  table  tests.  The  shake  table  test  results  show 
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the  stress  ratio  versus  the  number  of  cycles  leading  to  initial  liquefac¬ 
tion  (pore  pressure  in  the  specimen  equals  the  external  confining  pres¬ 
sure)  for  sands  at  four  different  relative  densities.  It  is  noted  that 
the  North  Island  data  plot  somewhere  above  the  =  68%  curve  for  the 
sand  in  the  low  cycle  range  and  above  the  54%  relative  density  line  for 
the  higher  cycle  range.  Thus,  although  the  North  Island  soil  appeared 
to  have  a  less  than  35%,  based  upon  quasi-static  and  dynamic  penetra¬ 
tion  tests,  it  exhibits  a  resistance  to  liquefaction  under  cyclic  triaxial 
conditions  comparable  to  that  of  a  denser  soil. 

This  increased  liquefaction  resistance  is  apparently  due  to  the 
distinct  structure  of  the  sensitive  silty  sandy  materials  in  layer  5. 

This  conclusion  is  further  supported  by  the  trend  in  the  test  data  in 
Figure  8.  It  is  noted  that  under  higher  stress  levels  (lower  number  of 
cycles  to  failure)  the  North  Island  soil  performs  like  one  with  a  higher 
relative  density  than  like  one  under  lower  stress  levels  (higher  number 
of  cycles  to  failure).  Possibly,  the  longer  time  to  failure  under  the 
lower  stress  levels  permits  more  disturbance  of  the  initial  soil  structure 
and,  hence,  causes  the  soil  to  perform  as  if  it  had  a  lower  relative 
density. 

Figure  9  portrays  some  of  the  results  of  a  compendium  of  soil  data 
compiled  from  a  number  of  commercial  testing  laboratories  (Ref  7).  The 
data  in  Figure  9,  for  both  uniformly  graded  (SP)  and  silty  sands  (SM), 
have  been  normalized  in  terms  of  D^.  (It  has  been  observed  that  by 
dividing  the  stress  ratio  causing  liquefaction  at  a  particular  stress 
ratio  by  relative  density,  reasonable  agreement  could  be  obtained  between 
the  liquefaction  resistance  of  similar  soils  at  different  densities,  as 
presented  in  Reference  14.  Also  shown  in  Figure  9  is  the  best  fit  curve 
(Ref  15)  for  data  from  a  number  of  laboratory  research  programs. 

It  is  noted  that  in  spite  of  any  differences  in  precision  between 
research  and  commercial  testing  activities,  the  mean  soil  strength 

curves  are  in  very  good  agreement/. 

/ 

Superimposed  upon  the  data  in  Figure  9  are  the  North  Island  soil 
test  data  from  Figure  7,  with  an  assumed  =  60%.  The  agreement  with 
the  mean  curves  indicates  that  the  North  Island  soil  had  a  resistance  to 
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liquefaction  roughly  equivalent  to  that  of  a  typical  soil  having  a 
relative  density  of  60%.  The  lower  values  estimated  for  this  soil 
(35%)  based  upon  the  penetration  readings  could  be  explained  in  terms  of 
the  sensitivity  of  the  soil  structure.  It  would  appear  that  the  soil 
structure  is  largely  lost  during  the  penetration  tests  in  advance  of  the 
sounding  device.  Thus,  the  penetration  tests,  both  the  split  spoon  and 
the  friction  cone,  do  not  recognize  the  inherent  resistance  to  cyclic 
liquefaction  provided  by  the  soil  structure  under  cyclic  triaxial  condi¬ 
tions.  Therefore,  the  most  critical  soil  stratum  in  the  profile  in 
Figure  1  might  better  be  characterized  as  that  having  the  strength  of  a 
sand  with  a  of  60%,  rather  than  the  35%  originally  assumed.  This,  in 
effect,  would  almost  double  its  liquefaction  resistance.  This  does  not 
necessarily  mean  that  the  soil  profile  in  question  does  not  still  repre¬ 
sent  a  liquefaction  hazard,  but  rather  that  the  threat  may  be  somewhat 
less  than  originally  concluded. 

It  must  be  kept  in  mind,  however,  that  should  sufficient  disturbance 
occur  to  destroy  the  in  situ  soil  structure,  the  soil  stratum  in  question 
could  perform  as  the  extremely  liquefiable  material  suggested  by  the 
penetration  results.  A  significant  question  that  comes  up  is  whether  or 
not  the  liquefaction  resistance  as  determined  by  the  laboratory  specimens 
is  most  valid  or  whether,  in  fact,  the  poorer  behavior  suggested  during 
the  penetration  tests  more  realistically  predicts  the  response  of  this 
area  under  an  actual  earthquake.  Unfortunately,  the  answer  to  this 
question  cannot  be  provided  at  this  time. 


SITE  EVALUATION  TECHNIQUES 


One  of  the  two  major  areas  being  addressed  in  this  program  is  that 
of  improved  reconnaissance  or  field  liquefaction  threat  evaluation 
techniques.  No  matter  how  sophisticated  and  reliable  analytical  proce¬ 
dures  may  become,  they  can  provide  accurate  response  predictions  only 


when  valid  input  data  can  be  obtained.  In  the  case  of  minor  structures, 
where  involved  analytical  programs  are  not  possible,  expedient  field 
evaluation  techniques  become  the  only  available  approach. 

The  liquefaction  potential  of  a  soil  is  directly  related  to  its 
volumetric  change  tendencies;  more  specifically  to  its  volumetric-strain/ 
shear-strain  coupling.  Traditionally,  the  volume  change  properties  have 
been  the  least  studied  from  the  standpoint  of  in  situ  measurements. 
Settlement  analysis,  predictions  of  swelling,  or  other  volume  change 
calculations  have  generally  been  based  upon  laboratory  testing  of  acquired 
samples.  Although  this  may  have  proven  adequate  in  the  past,  the  recent 
awareness  of  earthquake-induced  liquefaction  has  introduced  more  refined 
requirements.  Problems  associated  with  sample  disturbance,  laboratory 
simulation  of  in  situ  stress  state,  temperature,  chemical  and  biological 
environments,  and  soil  heterogeneity,  e.g.,  severely  impair  any  analysis 
based  upon  laboratory  testing. 

Thus,  although  in  situ  testing  is  attractive  in  theory,  direct 
determination  of  in  situ  volume  change  properties  applicable  to  earthquake- 
type  loading  is  not  yet  available.  Recourse  must  be  made  to  empirical 
correlations  between  volume  change  properties  and  some  form  of  expedient 
field  test.  Various  forms  of  penetration  tests  have  been  utilized  in 
this  regard,  particularly  the  standard  penetration  test  (ASTM  D  1586) 
and  the  friction  cone  penetration  test  (ASTM  D  3441-75T). 

Unfortunately,  the  penetration  resistance  of  a  soil  is  a  function 
not  only  of  its  volume  change  characteristics  but  also  of  its  strength, 
shear  stiffness,  and  other  deformational  characteristics.  Penetration 
resistance  can  be  influenced  by  factors  other  than  those  directly  influ¬ 
encing  liquefaction  potential.  Therefore,  it  may  be  necessary  to  measure 
several  types  of  response  before  accurate  penetration  test  correlations 
become  possible. 

One  promising  device  in  this  regard  is  a  piezometer  probe  (Ref  16) 
which  measures  pore  pressure  in  the  tip  of  a  soil  penetrating  cone.  If, 
during  cone  penetration,  positive  (increased)  pore  water  pressures  are 
generated,  then  the  effective  stresses  and,  hence,  the  strength  and 


resistance  to  penetration  arc  reduced.  Alternat  tvely ,  il  negative 
l reduced)  pore  water  pressures  occur,  then  the  soil  structure  is  dilating, 
effective  stresses  are  increased  and  so  is  penetration  resistance.  I'lie 
measured  incremental  changes  in  pore  pressure  during  penetration  ot 
saturated  soils  are  directly  related  to  volume  change  tendency.  They 
are  also  a  function  of  soil  permeability  and  rate  ot  penetration. 

The  rate  at  which  pore  water  pressures  reach  o<pii librium  tol lowing 
cessation  ot  penetration  are  direct  functions  ot  permeability.  I'lms ,  a 
device  exists  with  the  potential  lor  correlation  with  hotli  volume  change 
characteristics  and  permeability.  These  two  tactors  (along  with  nature 
of  loading)  are  the  major  determinants  which  control  the  occuiieuce  and 
severity  of  soil  ! i quo  tact  ion. 

A  schematic  ot  the  piezometric  probe  (Ket  lb)  developed  hv  Wissa  is 
shown  in  Figure  10.  The  rapid  response  is  made  possible  by  means  ot  a 
Ingh  air  entry,  stainless  steel,  porous  cone  tip  hydraulically  connected 
to  a  flush  diaphragm  pressure  transducer.  A  four-conductor  shielded 
cable  encased  in  polyethylene  tubing  communicates  with  the  surtace 
through  standard  "A"  drill  rods. 

Figures  11  and  12  show  typical  piezometric  output  readings.  Figure 
II  shows  differences  in  pore  pressure  generation  and  dissipation  rates 
among  tour  different  types  ot  soil. 

It  is  noted  that  loose  soils  l  and  2  lend  to  compact  during  yielding 
caused  by  cone  penetration  and,  hence,  generate  increased  pote  watei 
pressures.  These  induced  pore  water  pressures  then  dissipate  a(  two 
different  rates:  the  sand  in  a  matter  ot  a  few  minutes  and  the  clayey 
silt  over  a  period  ol  almost  2  hours. 

Denser  soils  i  and  ■)  dilate  dining  penet  rometer- induced  yielding, 
causing  reduction  in  pore  water  pressure*  below  that  ot  the  datum. 

Again,  stabilizing  ol  pore  water  pressure  is  a  function  ot  permeability, 
with  silt  requiring  about  10  times  longer  tor  pressure  equalization  t ban 
dense  sand  requires. 

A  comparison  ot  the  behavior  ot  similar  soil  types  at  different 
densities  is  shown  in  Figure  12.  As  hypothesized,  soils  in  the  dense 
state  dilate,  whereas  those  m  the  loose  state  compact.  The  time  ot 
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pore  pressure  equilization  is  shown  to  he  both  a  function  of  soil  type 
(grain  size)  and  state  of  initial  compactness;  (i.e.,  loose  or  dense). 

Obviously,  substantial  differences  exist  between  stresses  induced 
hv  seismic  events  such  as  earthquakes  and  stresses  that  occur  during 
cone  penetrations.  Nevertheless,  since  responses  in  both  cases  are  a 
function  of  soil  volume  change  tendency  and  permeability,  useful  correla¬ 
tions  may  be  possible. 

An  entirely  different  approach  to  liquefaction  prediction  is  one 
based  upon  soil  resistivity  measurements  (Ref  17).  The  liquefaction 
potential  of  a  saturated  soil  is  largely  influenced  by  the  state  of 
compactness  (or  relative  soil  density).  However,  many  additional  factors 
can  have  an  influence  here.  One  term  used  in  conjunction  with  the 
strength  of  cohesionless  soils  is  "fabric."  This  refers  to  the  structure, 
or  the  way  in  which  the  individual  grains  interlock  or  fit  together. 

The  fabric  of  a  soil  is  largely  controlled  by  the  way  in  which  the  soil 
was  originally  formed  and  its  subsequent  stress  or  deformationa 1  history. 
This  factor  has  been  investigated  recently  by  considering  different 
types  of  sample  preparation  and  conditioning  (Ref  18-21).  Substantial 
results  now  exist  to  support  the  fact  that  soil  fabric  exerts  a  measurable 
effect  on  soil  liquefaction  potential,  at  least  in  the  laboratory. 

Thus,  if  it  were  possible  to  measure  soil  density  and  fabric  in  situ,  it 
is  conceivable  that  greatly  improved  liquefaction  prediction  capability 
would  be  available. 

A  device  being  applied  by  Arulanandan  (Ref  17)  appears  to  have 
considerable  potential  in  this  regard.  This  approach,  by  measuring 
conductivity  in  two  directions,  attempts  to  characterize  the  soil  struc¬ 
ture  to  account  for  shape  and  spatial  arrangement  of  particles.  Thus, 
void  ratio  and  a  measure  of  structural  anisotropy  can  be  deduced  which 
can  be  related  to  cyclic  loading  behavior. 

The  approach  utilizes  a  parameter  called  the  formation  factor,  K, 
defined  as  the  ratio  of  conductivity  of  the  electrolyte  (pore  water)  to 
the  conductivity  of  the  sand  saturated  by  the  electrolyte.  Considerable 


il.it.  i  exists  (Hi't  22  ami  2.1)  supporting  the  tol  lowing  relationship  between 
porosity,  n,  and  F: 


Other  work  (Kef  24)  has  extended  this  relationship  to  un  hide  a  part  tele 
shape  factor,  x,  which  is  a  measure  ot  the  sphericity  of  soil  grains. 

For  example,  x  =  2  for  perfect  spheres  and  decreases  as  the  axial  ratio 
ot  the  spheroid  increases.  This  relationship  is  of  the  form 


F  =  (x  *  D  -  n 
xn 

The  problem  with  such  relationships  is  that  a  correlation  foi  the 
particular  sand  had  to  be  established  tiist  (generally  in  the  laboratory), 
and  this  failed  to  preserve  the  initial  in  situ  fabric.  More  recent 
developments  at  the  University  of  California  at  Davis  (Ket  17)  are 
designed  to  bypass  this  problem.  This  work  hypothesizes  that  the  tesis- 
tivity  path  involved  in  measurement  of  F  is  a  function  ot  particle 
orientation.  Thus,  by  measuring  an  anisotropic  index,  A  (defined  as  the 
square  root  ol  the  ratio  of  vertically  measured  to  horizontally  measured 
formation  factors)  in  situ,  a  measure  of  the  natural  soil  fabric  can  be 
measured,  prior  to  any  sampling  disturbance.  There  are  still  problems 
to  be  addressed  such  as  the  effects  ol  any  interparticle  cementations 
and  assessing  of  the  in  situ  lateral  pressures,  Nevertheless,  this 
approach  does  show  promise  of  providing  valuable  in  situ  definitive  soil 
data . 

Various  other  devices  and  procedures  are  available  tor  determining 
m  situ  volume  change  tendencies  and,  hence,  liquefaction  potential. 

These  include  borehole  permeability,  soil  pressuremet er,  cvclii  screw 
plate,  and  other  tests;  however,  these  test  methods  are  not  considered 
suitable  for  further  discussion  herein.  The  borehole  permeability  test 


concerns  itsoll  primarily  with  long-term  consolidation  effects.  The 
other  two  tests,  even  it  combined  with  pore  pressure  readings,  would  be 
difficult  to  interpret  from  a  liquefaction  evaluation  viewpoint.  Other 
more  esoteric  techniques  are  available,  such  as  the  water  cannon  (Rel 
25)  developed  at  the  Technical  University  at  Zurich  or  the  cylindrical 
in  situ  (CIST)  test  (Ref  26),  but  such  tests  are  not  considered  compati¬ 
ble  with  the  aims  of  this  research  program. 


COMPUTER  MODELS 

Advanced  Analytical  Models 

Significant  research  is  in  progress  to  model  the  generation  and 
dissipation  of  pore  pressure  with  seismic  excitation.  Present  implemented 
technology  is  essentially  limited  to  total  stress  analysis  (i.e.,  ignoring 
the  two-phase,  soil-pore  fluid  interaction):  however,  current  research 
is  studying  effective  stress  techniques.  The  transition  from  both  total 
stress  analysis  and  effective  stress  analysis  to  a  tool  that  goetechnical 
engineers  can  use  to  solve  practical  design  problems  is  not  imminent. 

The  concept,  however,  represents  a  critical  advance  which  could  be  very 
s ign i f icant  . 

Tot  a  1  Stress  Analysis 

Several  techniques  are  available  for  total  stress  analysis.  SHAKE 
(Ref  27)  and  APOLLO  (Ref  28)  (also  GADFLEA,  Ref  29)  have  been  used 
together.  SHAKE  generates  total  stress  histories  that  are  converted  to 
equivalent  stress  cycles.  These  data  and  the  soil  strength  are  used  to 
determine  the  undrained  number  of  cycles  to  liquefaction,  which  is  input 
to  APOLLO.  A  one-dimensional  code,  APOLLO,  and  a  two-dimensional  code, 
GADFLEA,  compute  the  uncoupled  pore  pressure  buildup  and  dissipation. 

No  provision  is  included  to  account  for  variation  of  soil  modulus  with 
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increasing  pore  pressure  and  reduced  effective  confinement.  SHAKE  uses 
equivalent  linear  soil  properties,  which  are  adjusted  for  strain  ampli¬ 
tude.  An  example  is  illustrated  below. 

Narasinthan  (Ref  30)  uses  a  seepage  program  -  TRUST  -  to  analyze  the 
dissipation  phase  of  pore  pressure  production.  The  initial  pressures 
are  estimated  by  undrained  total  stress  analysis. 

Effective  Stress  Analysis 

Finn  et  al.  (Ref  31)  have  developed  DESRA  which  is  a  one-dimensional, 
effective-stress,  soil  model  that  computes  pore  pressure.  This  is  a 
nonlinear  procedure  in  which  the  equations  of  motion  are  integrated 
using  the  Newmark  method.  The  initial  loading  is  given  by  a  hyperbolic 
stress-strain  relation.  Reloading,  described  by  the  Masing  criterion, 
produces  hysteresis  loops.  The  upper  curve  of  the  hysteresis  loop  can 
be  obtained  from  the  loading  curve  by  translating  the  origin  and  increas¬ 
ing  the  scale  vertically  and  horizontally  by  a  factor  of  two.  When  used 
in  the  effective  stress  mode,  DESRA  produces  increases  in  pore  pressure 
with  each  cycle,  using  the  parameters  which  describe  the  volume  change 
and  rebound  characteristics  of  the  soil.  The  reduction  in  effective 
stress  results  in  reduction  of  shear  modulus.  It  includes  the  effects 
of  both  shear  strain  and  progressive  increase  in  pore-water  pressure  on 
soil  properties.  However,  it  is  limited  to  one-dimensional  cases  as 
well  as  equivalent  cyclic  input  (not  actual  time  histories),  and  it  does 
not  treat  dissipation. 

Ghaboussi  and  Dikmen  (Ref  32)  have  a  program  of  similar  capabilities. 
This  program,  although  one-dimensional,  uses  actual  time  history  input 
as  the  loading  function.  The  program  utilizes  a  material  model  similar 
to  that  of  Ishihara  (Ref  33  and  34),  modified  by  an  elliptic  stress 
loading  path.  The  program  treats  only  effective  stress  generation  and 
not  dissipation. 
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Baladi  and  Rohani  (Ref  35)  have  modified  the  "cap  model"  to  pore 
pressure  generation.  This  model  is  still  in  the  research  stage  and  has 
not  been  implemented. 


Total  Stress  Analysis  Example 

Seed  et  al.  (Ref  28)  have  investigated  the  distribution  of  hydro¬ 
static  pore  pressure  u  in  the  soil  by  use  of  the  one-dimensional  equation 


where  C  =  coefficient  of  consolidation  of  the  soil 

v 

z  =  depth  within  soil  (the  vertical  coordinate) 

3u^/3t  =  rate  of  pore  pressure  generation  caused  by  earthquake 

This  is  the  diffusion  equation  used  in  Terzaghi's  classical  consolidation 
theory,  with  a  pressure-generating  term  added.  The  solution  of  this 
equation  is  accomplished  by  the  finite-difference  technique  using  incre¬ 
mental  time  steps.  The  pore  pressure  generation  is  estimated  by  Figure  13 
as  a  function  of  the  number  of  cycles  to  cause  liquefaction. 

The  coefficient  of  consolidation  C  ,  which  is  defined  in  terms  of 

v’ 

the  coefficient  of  volume  compressibility  m  and  the  coefficient  of 
permeability  k,  may  be  estimated  by  means  of  Figures  14  and  15. 


where  y^  =  the  weight  of  water 

The  rise  in  the  water  table  All,  during  time  increment  At,  is  given  by: 


AH 


n 

t' 


where  n  =  the  effective  porosity 

This  procedure  has  been  automated  in  the  computer  program  APOLLO 

and  may  be  used  with  the  analysis  in  the  computer  program  SHAKE.  SHAKE 

is  used  to  produce  the  equivalent  uniform  cyclic  shear  stress,  T  ,  and 

eq 

the  equivalent  number  of  uniform  stress  cycles,  n  ,  for  various  depths 

eq 

ot  soil.  From  strength  data  the  number  of  cycles  to  cause  liquefaction 
at  each  depth  is  determined.  Using  this  information,  program  APOLLO 
solves  the  pore  pressure  generation-dissipation  equation. 

1’he  pore  pressure  generation  function  is  based  on  undrained  test 
data.  This  application  is  deemed  sufficiently  accurate  when  small  time 
steps  are  used  to  properly  account  for  drainage.  The  elastic  response 
analysis  used  to  determine  the  number  of  cycles  to  liquefaction  can  be 
made  in  order  to  consider  the  isolation  effects  of  subsurface  liquefaction 
on  near  surface  shaking  and  the  reduction  in  pore  pressure  generation 
when  iteration  techniques  are  used. 

A  typical  example  from  Seed  et  al.  (Ref  28)  from  the  Niigata  earth¬ 
quake  of  I9b4  is  shown  in  Figures  16  and  17.  The  computed  variations  ot 
pore  water  pressure  with  time  are  given.  Figure  17  shows  the  buildup  of 
pore  pressures.  It  may  be  seen  that  the  sand  layer  at  a  depth  of  15 
feet  liquefies  after  about  21  seconds  of  shaking;  liquefaction  extends 
to  depths  of  20,  30,  and  40  feet  after  about  23,  32,  and  40  seconds  of 
shaking.  Although  the  layers  above  a  15-foot  depth  continue  to  increase 
in  pore  pressure  as  the  shaking  progresses,  the  rate  of  increase  is  very 
low  after  the  15-foot  level  liquefies.  Reference  2  notes  that  when  the 
pore  pressure  ratio  in  the  top  foot  of  soil  reaches  bl!%,  the  ground  will 
become  soft,  and  a  man  will  sink.  This  occurred  after  about  8.5  minutes 
in  Niigata,  determined  in  the  analysis.  The  pore  pressure  ratio  at  the 
ground  surface  begins  to  decrease  after  about  20  minutes  but  will  not 
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support  a  man  until  about  40  to  50  minutes  after  the  earthquake.  The 
results  of  the  computer  analysis  are  in  general  agreement  with  observed 
reports . 

If  the  water  table  were  located  at  a  depth  of  15  feet,  no  significant 
pore  pressure  increases  would  occur  in  the  upper  10  feet  of  soil,  even 
though  the  soil  is  liquefied  between  15  and  40  feet.  Thus,  in  this 
situation  the  hearing  capacity  of  small  shallow  footings  near  the  surface 
might  well  be  essentially  unaffected  by  the  dissipation  of  pore  water 
pressures  in  the  liquefied  zone. 

Program  APOLLO  has  been  expanded  into  a  two-dimensional  computer 
program  called  GADFLEA.  The  approach  is  very  similar  to  the  one- 
dimensional  analysis,  requiring  as  input  information  the  number  of 
cycles  causing  liquefaction  by  soil  element.  The  number  of  cycles 
causing  liquefaction  is  a  function  of  the  applied  shear  stress  loading 
and  soil  confinement.  These  may  be  determined  from  a  conventional 
two-dimensional  elastic  or  inelastic  finite  element  analysis.  Using  the 
input  data  program,  GADFLEA  computes  the  two-dimensional  pore  pressure 
generation  and  dissipation  from  the  earthquake. 

Effective  Stress  Analysis  Example 

Ghaboussi  and  Dikmen  (.Ref  32),  using  their  effective  stress  model, 
also  studied  the  Niigata  site.  Ghaboussi  uses  two  distinct  material 
behavior  models  to  cover  the  range  from  the  initial  in  situ  stress 
condition  to  liquefaction.  The  first  model  represents  the  soil  behavior 
from  the  in  situ  stress  condition  up  to  the  initial  liquefaction.  A 
second  material  model  is  used  to  treat  post-liquefaction.  The  initial 
liquefaction  is  to  be  considered  a  condition  of  near  failure  (.not  100*, 
pore  pressure  confining  stress  ratio). 

Preliquefaction.  The  shear  stress-strain  relationship  under  mono¬ 
tonic  loading  is  represented  by  the  following  equation  and  shown  in 
Figure  18a. 


17 


s 

max 


in  which  q  is  the  shear  stress,  p'  is  the  effective  pressure,  and  y  is 
the  shear  strain.  The  factors  S  and  G  are  defined  in  the  figure. 

Tlie  unloading  takes  place  linearly  with  the  slope  until  previous 
maximum  or  minimum  value  of  q/p '  is  reached,  whereupon  the  above  stress- 
strain  relationship  becomes  valid  again.  The  relationship  is  active  up 
to  the  onset  of  initial  liquefaction.  Experimental  evidence  (Ref  33) 
suggests  that  in  the  p'-q  plane  the  shear  yield  loci  (lines  of  equal 
shear  strain)  take  the  form  of  straight  lines  radiating  from  the  origin 
(Figure  18b).  At  higher  values  of  effective  pressure,  the  yield  loci 
approach  a  state  parallel  to  the  effective  stress  (p ' )  axis.  Under 
monotonic  shear  stress  increase,  failure  occurs  on  a  line  corresponding 
to  an  asymptote  of  the  shear  stress  relationship  indicated  earlier. 

This  critical  state  line  fj  is  given  by 

fj  =  q  -  p'  tan  ip  =  0 

where  tan  i }»  =  the  frictional  strength  of  the  soil 

The  stress  path  to  failure  i,  has  been  approximated  by  a  quarter  of  an 
ellipse  given  by  the  following  equation. 

f2  =  (p'  '  Pp2  +  «T  ’  (Pq  '  PfT  =  0 


where  the  subscript  f  is  used  to  designate  conditions  of  failure.  The 
material  parameter  \  is  the  ratio  of  the  major  and  minor  axes  of  the 
ellipse,  as  given  by  the  following  relation. 


$sp?SP3-<£' iit: 


A 


0 


The  stress  path  f.,  is  completely  defined  by  the  two  material  parameters 
A  and  0  and  the  in  situ  effective  pressure  p^  and  can  be  written  in 
terms  of  the  following  parameters. 


The  effective  pressure  p'  decreases  as  a  result  of  increases  in 
pore  pressure.  In  unloading,  the  effective  pressure  remains  constant 
until  the  stress  path  crosses  previous  maximum  or  minimum  value  of  a/p' , 
at  which  point  a  new  value  of  p^  is  determined,  using  the  f.,  function. 
The  parameter  A  is  related  to  relative  density  (Figure  19).  Figure  19 
is  based  on  fairly  uniform  sands  with  rounded  particles.  Variation  in 
shape  has  some  effect  on  A. 


Post  Initial  Liquefaction.  The  material  behavior  of  the  soil 
changes  abruptly  after  initial  liquefaction.  The  stress  path  is  defined, 
based  on  very  limited  soil  data,  but  further  development  is  required  in 
this  area  of  investigation.  After  the  earthquake  motion  has  ceased, 
excess  pore  pressures  exist  in  the  soil  layers.  These  dissipate  with 
time  and  require  analysis  using  a  dissipation  program. 

Niigata  Case  Study.  The  Niigata  site  discussed  earlier  was  studied 
by  Ghaboussi  (Ref  32).  Figures  20,  21,  and  22  give  results  of  that 
analysis.  The  effective  stress  analysis  shows  that  the  increase  in 
effective  pressures  with  depth  have  a  significant  influence  on  the 
maximum  shear  stresses,  equivalent  shear  moduli,  and  hysteretic  damping. 
Although  modelling  differences  do  exist,  the  qualitative  behavior  is 
very  similar  to  that  of  the  total  stress  analysis  in  predicting  the  zone 
that  will  undergo  liquefaction.  However,  an  effective  stress  analysis 
allows  for  strain  softening  not  possible  with  a  total  stress  analysis. 
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This  may  be  significant  when  isolation  from  shaking  of  upper  layers 
occurs  as  a  result  of  weakened  (liquefied)  deeper  layers.  Though  much 
additional  work  is  required  on  this  model,  it  appears  to  have  great 
potential.  It  is  presently  limited  to  one-dimensional  preliquefaction 
ana  lysis. 


SUMMARY  OF  PAST  EFFORT 

Technology  for  defining  the  soil  liquefaction  phenomenon  is  devel¬ 
oping  through  many  different  approaches.  This  is  in  response  to  the 
knowledge  that  soils  are  a  critical  constituent  of  most  civil  engineering 
structures.  As  indicated  in  Reference  36,  poor  soil  behavior  during 
past  earthquakes  has  led  to  structure  failures.  The  following  are 
examp les . 

o  Relatively  moderate  ground  shaking  in  the  1964  Niigata, 

Japan,  earthquake  caused  extensive  ground  failure  and  loss 
of  support  for  buildings  due  to  soil  liquefaction,  which 
resulted  in  damage  approaching  $1  billion. 

o  The  catastrophic  landslides  and  other  soil  failures  in  the 
1964  Alaska  earthquake  destroyed  much  of  the  developed  resi¬ 
dential  and  commercial  property  in  several  cities  as  well 
as  many  highway  bridges. 

o  During  the  1906  San  Francisco  earthquake,  ground  cracking 
and  sliding  caused  rupture  of  many  buried  water  lines,  thus 
immobilizing  fire  fighting  operations  resulting  in  uncon¬ 
trolled  fire  sweeping  through  the  city. 

o  lu  1972  an  earthquake  in  Peru  produced  an  enormous  avalanche 
and  landslide  which  rushed  down  the  mountain  and  buried  entire 
villages  in  the  valley  below. 
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In  addition,  two  earth  dams  nearly  tailed  in  tlie  1971  San  Fernando 
earthquake.  This  failure  would  have  inundated  a  residential  area  of 
80,000  people  and  could  have  been  the  greatest  single  natural  disaster 
in  the  history  of  our  country. 

Although  all  ot  these  tragedies  cannot  be  attributed  solely  to  the 
liquefaction  process,  it  is  widely  accepted  that  soil  liquefaction  is 
the  major  earthquake-associated  threat  confronting  geotechnical  engineer¬ 
ing.  In  spite  of  the  increased  efforts  expended  for  solution  of  this 
problem,  it  was  the  consensus  reached  at  a  1977  workshop  attended  by  72 
internationally  recognized  experts  (Ref  36)  that  available  geotechnical 
earthquake  engineering  techniques  and  design  procedures  are  either 
questionable  or  inadequate.  This  evaluation  is  particularly  pertinent 
to  the  watcrtiont  area.  Marine  earthworks,  including  such  waterfront 
structures  as  tied  bulkheads,  cofferdams,  rockfill  jetties,  and  partially 
relieved  quay  walls,  are  commonly  situated  in  relatively  precarious 
regions.  Often,  the  soil  fill  associated  with  such  structures  is  poorly 
compacted.  Evaluation  techniques  are  necessary  for  assessing  (1)  the 
integrity  of  soil  foundations  for  waterfront  structures  and  (2)  in 
particular,  the  susceptibility  of  the  fill  materials  to  loss  of  stability 
during  seismic  loading. 

The  research  program  described  herein  is  concerned  with  technology 
transfer,  as  well  as  with  individual  or  unique  developments.  As  such, 
it  concentrates  primarily  upon  two  areas:  (1)  expedient  and  reliable 
site  evaluation  techniques  and  (2)  valid  analytical  response  prediction 
procedures . 

To  date,  an  appraisal  of  the  best  available  analytical  computer 
codes  for  evaluating  liquefaction  has  been  conducted.  The  Navy  is 
offering  support  to  the  most  promising  ideas. 

A  report  has  been  published  (Ref  4)  which  offers  preliminary  guidance 
to  the  Naval  Shore  Establishment  until  more  adequate  liquefaction  tech¬ 
nology  is  available.  This  report  summarizes  the  current  state-of-the-art 
regarding  placement  of  structures  in  liquefiable  regions.  Field  studies 
were  continued  at  NAS  North  Island,  which  had  previously  (Ref  2)  been 
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identified  as  representative  of  a  liquefaction  hazard.  Undisturbed 
samples  of  waterfront-unique,  potentially  liquefiable  soils  were  obtained 
and  tested  in  the  laboratory  under  simulated  earthquake  conditions 
(Ref  13).  The  information  secured  indicated  -  at  least  for  the  one 
highly  structured  coastal  soil  tested  to  date  -  that  liquefaction  resis¬ 
tance  appears  to  be  somewhat  greater  than  that  based  on  correlations 
derived  from  simple  field  tests. 


FUTURE  WORK 

Work  for  the  upcoming  year  will  be  directed  along  three  different 
avenues  of  study,  with  roughly  equal  emphasis:  (1)  improvement  of 
analytical  models,  (2)  investigation  of  innovative  site  reconnaissance 
techniques,  and  (3)  more  detailed  studies  of  soil  parameters  with  parti¬ 
cular  attention  to  the  Navy's  unique  soils,  using  CEL's  newly  acquired 
liquefaction  testing  capabilities. 

With  regard  to  the  latter  area,  only  one  of  the  problem  soils  at 
NAS  North  Island  (a  very  sensitive  coastal  silt  deposit)  has  been  sub¬ 
jected  to  cyclic  laboratory  testing  to  date  and  appears  to  pose  a  major 
liquefaction  threat.  Several  other  unusual  soils  exist  at  this  facility 
One  soil  of  interest  is  a  hydraulic  fill  sand  containing  an  unusually 
high  percentage  of  mica  flakes.  How  this  type  of  dredged  fill  material 
will  behave  under  earthquake  shaking  must  be  determined.  Some  forms  of 
parameter  studies  of  these  coastal  marine  deposits  must  be  carried  out. 
For  example,  it  has  been  observed  that  when  a  substantial  portion  of  any 
granular  soil  passes  the  #200  mesh  sieve  the  penetration  resistance  is 
drastically  reduced  in  both  static  and  dynamic  tests.  The  significance 
of  this  phenomenon  should  be  determined  for  actual  dynamic  or  cyclic 
behavior . 

It  is  intended  that  the  most  promising  devices  such  as  piezometric 
cones  or  directional  resistivity  measuring  instruments  be  evaluated  at 
the  North  Island  site. 
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With  regard  to  analytical  developments,  it  is  intended  that  the 
most  promising  available  soil  models  will  be  incorporated  into  CEL's 
library  of  finite  element  codes,  which  would  then  be  exercised  on  parti¬ 
cular  solutions.  Eventually  such  a  code  would  be  combined  with  a  ground- 
water  flow  code  (such  as  TRUST)  for  complete  definition  of  soil-structure 
interaction  under  dynamic  and  hydrodynamic  conditions. 
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Figure  A.  Correlations  between  vertical  effective  stress,  relative 
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Relationship  between  shear  stress/confining  stress  ratio 
lOhv^av^o^  ^nd  relative  density  for  known  cases  of  lique 
faction  and  nonliquefaction  (Ref  10). 
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Figure  9.  Normalized  stress  versus  number  of  cycles  to  initial  liquefaction  (Ref  7). 
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Time-rate  of  dissipation  of  excess 
pore  pressure  generated  by  the  probe 
in  dense/stiff  and  loose/soft  soils 
(Ref  16). 


Cyclic  Ratio 


Figure  13.  Rate  of  pore  water  pressure  buildup  in 
cyclic  simple  shear  tests  (Ref  10). 


Peak  Pore  Pressure  Ratio 


Figure  14.  Theoretical  relationships  between  compres¬ 
sibility  of  sands  and  pore  pressure  buildup 
(Ref  28). 
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Pore  Pressure  Ratio,  u/o' 


(a)  During  earthquake  shaking. 
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(b)  In  8-minute  period  following  earthquake. 

Figure  17.  Computed  development  and  variation  of  pore 
water  pressures  for  soil  profile  shown  in 
Figure  16  (Ref  28) . 
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(c)  In  60-minute  period  following  earthquake. 
Figure  17.  Continued. 
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(a)  Normalized  stress- 
strain  relationship. 


(b)  Stress  path. 


Figure  18.  Ishihara  material  model. 
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Figure  19.  Relationship  between  material  parameter  X  and  relative 
density. 
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Figure  20. 


Variation  of  pore  pressure  with  depth  (Ref  32) 
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Acceleration 
Peak  acceleration 

Coefficient  of  consolidation  of  the  soil 
Relative  density  of  the  soil 

Formation  factor,  used  to  characterize  a  soil  structure 

Constants  in  soil  model  from  Reference  32 

Unloading  shear  modulus 

Thickness  of  soil  layer 

Soil  permeability 

Soil  coefficient  of  volume  compressibility 
Number  of  cycles  to  liquefaction 
Soil  porosity,  including  non-interconnected  voids 
Effective  soil  porosity 

Equivalent  number  of  stress  cycles,  used  in  conjunction 
with  X 

eq 

Effective  confining  pressure 
Effective  pressure  at  failure 
Initial  effective  pressure 
Shear  stress 

Maximum  possible  shear  stress  to  effective  stress  ratio 

Symbol  for  a  silty  sand 

Symbol  for  a  poorly  graded  sand 

Pore  water  pressure 

Shape  factor  used  to  characterize  particle  shape 
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LIST  OF  SYMBOLS 
(Continued) 


Depth  coordinate  distance 


Increment  of  pore  water  pressure 


Lateral  confining  pressure 


Parameter  used  for  soil  model  in  Reference  32 


Ratio  of  maximum  shear  stress  to  confining  pressure 


Initial  effective  volumetric  stress 


Average  horizontal  shear  stress 


Time  derivative  of  rate  of  pore  water  pressure  generation 
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SOU  I  II WI-'S I  RSCII  INSI  King.  Snn  Antonio.  I  X.  R  IVHutl,  Sun  Antonio  I'X 


SI  ANFORD  UNIVI  RSI  TY  Engt  I  ib.  Stanford  CA;  Stanford  <  A  (GcncJ 
SIAM  UNIV.  Ol  NEW  YORK  Buffalo.  NY 

TEXAS  AAM  UNIVI  RSI  I  V  College  Siauou  I  X  (l  l  Dept  Herbich);  W  II  I  edbetlei  College  Siafion.  IX 
UNIVI  RSIIV  Ol  UAI  IEORNIA  1IERKET  IV.  CA((I  Dl  PI .  GERWICK)  HI  RKI  I  I  Y.CAiCI  HI  I’l . 
MIIVHI  l  I  I.  Berkeley  UA  (II  Hresler);  Berkeley  C’A  ilVpl  ol  Naval  Arch.l.  Beikeley  UA  (R  Williamson): 
DAVIS.  CAtCE  DI  PT.  I  AVT  OR);  LIVERMORE.  CA  (l  AWRINCI  LIVERMORE  I  AH.  TOKAR/I.  I  a  folia 
UA  f  Vcq  IVpl.  I  ib  C-07.VA);  M  Ounean,  Berkeley  UA;  SAN  DIEGO.  UA.  I  A  JOI  I  A.  UA  rS|  RiK  kll 
UNIV  I  RSI  I  V  Ol  Ol  I  AWARE  Newark,  01  (IVpl  of  Civil  I- nitincei mgt.  Chesson) 

UNIVI  KSI  I  \  OE  HAWAII  HONOLULU.  Ill  (SCIENCE  ANDTT  CH  OIV.I;  Honolulu  HI  ilfi  S;il.mll 
UNIVI  RSIIV  Ol  II  I  INOIS  Mel/  Rel  Rm.  Urbana  II  .  URBANA.  It.  (DAVISSON I;  URBAN  A.  II  il  IHR  VRV  I. 

URBANA.  II  i N I  VV  MARK);  Urbana  II  (UE  IVpl.  W  Gamble) 

UNIVI  RSIIV  Ol  M  ASS  ACHUSI  I  1 S  (Heroneimis),  Amherst  MAUI-  Oepl 

UNIVI  RSIIV  Ol  MICHIOAN  Ann  Arbor  Ml  (RichnrO 

UNIVI  RSIIV  Ol  NEBRASKA  I  INCOI  N  Lincoln,  Nl  (Ross  Ice  Shell  Proj  l 

UNIVI  RSI  I  V  Ol  Nl  VV  Ml  XICO  .1  Nielson  I- njjr  Malls  A  Civil  Sy  s  On  .  Albuquerque  NM 

UNIVI  RSIIV  Ol  TEXAS  Inst  Marine  Set  (l  ibrary).  Port  Arkansas  TX 

UNIVTRSI IV  Ol  TEXAS  AT  AUSTIN  AUSTIN.  TX  (  THOMPSON).  Austin  VX  fR  Olson);  Austin.  TX  (Breen) 
UNIV  I  RSII  V  1)1  VV  ASHINGTON  IVpl  of  Cisil  I  ngr  (Di  Malioekl.  Seattle  WA.  SEATII  E.  WA  (MERCHANT); 
SI  VTTI  E,  WA  (OCEAN  ENG  RSCH  LAB.  GRAY);  Seattle  VV  A  (I-  1  inger).  Seattle.  VV  A  I  lanspoilalion. 
Construction  «Y  (ieoin.  On 

VIRC.IN1A  INST  Ol  MARINI  SCI  Gloucester  Point  VA  (Library) 

Al  I  Rl  O  A  VII  .V  ASSOC  Honolulu  HI 

AMEIEK  ( tffshore  Res  A  Engl  Oiv 

Al’PI  110  IT  CH  UOUNCH  R  Scholl.  Palo  Alto  CA 

ARVIO  GRANT  Ol  YMPIA.  WA 

All  ANTIC  RICH  Eli  I  DUO.  0A1  I  AS.  TX  (SMI  I  ID 

AUSIRAI  IA  IVpl  I’W  (A.  Hicks).  Melbourne 

III  CHIT  I  CORP.  SAN  TRANCISCO.  CA  (PHI-UPS) 

BELGIUM  HAEUON.  N.V..  Gem 

III  ITII  MU  M  SIT- 1  I  CO  Oismuke.  Bethelehem.  PA 

IIOUW  KAMI’  INC  Berkeley 

BROWN  A  CAI.DWM  I  I  M  Saunders  Walnut  Creek,  CA 
BROWN  A  ROOT  Houston  TX  (0.  Ward) 

CANADA  Mem  Unis  Newfoundland  (Chari).  St  Johns.  Surveyor.  Nenninger  A  Chcnevert  Inc  .  Montreal;  Trans  Mm 
Oil  Pipe  I  one  Coip.  Vancouver.  BU  Canada;  Warnoek  Hersev  Prof  Srv  I  id.  I  a  Sale,  fjuebec 
a  BR  AUN  CO  On  Bouehct,  Murray  Hill.  N.I 

CHEVRON  OH  Ml- 10  Rl  SI  ARCH  CO.  I  A  HABRA.  CA  (BROOKS) 

CONCRETE  TECHNOl  ( KiY  CORP  TACOMA.  WA  (ANDERSON) 

CONRAD  ASSOC  Van  Nuys  CA  (  A  I  nisoni) 

ORAVO  CORP  Pittsburgh  I'A  (Wright) 

NORWAY  01  I  NORSK  I'  VI  Rl  TAS  (library).  Oslo 

EVAI  UATION  ASSOC.  INC  KING  OK  PRUSSIA.  PA  (IT  Oi  l  K.l 

I  ORD.  BACON  A  DAVIS.  INC  New  York  (Library) 

l-RANCE  Ol  .  Outertie.  Boulogne;  I  Pliskin,  Paris;  P.  Jensen.  Boulogne 

GKOTK.UHNICAI .  KNGINKKRS  INC  Winchester.  MA  (Paulding) 

GI.IDDEN  CO  STRONGSVll  l.l  .  OH  (RSCH  I  IH) 

GRUMMAN  AEROSPACE  CORIV  Bclhpage  NY  (Tech  Info  Ctrl 
HAI  I  V  A  Al  ORICH.  INC  Cambridge  MA  (Aldrich.  Jr.) 

HONEYWEl  I  .  INC  Minneapolis  MN  (Residential  Engt  I  ib.) 

HUGHES  AIRCRAFT  Culver  City  CA  (  lech.  Doe.  Cii ) 

ITALY  M,  Cninrni,  Milan;  Sergio  Tatloni  Milano;  Torino  (E.  I  evi) 

MAKAI  OCEAN  ENGRNG  INC  Kailua.  III 
JAMES  CO.  R  Girdley.  Orlando  FL 

I  .AMONT-DOHERTY  GEOLOGICAL  OHSERV  Palisades  NY  (McCoy);  Palisades  NY  (Sclwvn) 

LOCKHEED  MISSILES  A  SPACE  CO  INC.  Mgr  Naval  Arch  A  Mar  Eng  Sunnyvale.  CA;  Sunnyvale  UA 
(Rynewicr);  Sunnyvale,  UA  (K.L.  Kerr) 

LOCKHEED  OCEAN  LABORATORY  San  OicgoCA  (I  Simpson) 

MARATHON  OIL  CO  Houston  TX 


— — - -  ^  — ***jf-- — ■  -  - -  '"poi1  -‘w.'-  wi-””  mjfi  «. 


Men  II  I  AND  I  NtilNI  IKS  INI'  Houston  IX  (B  McClelland  I 
Ml  IXJNNIL  AIRCRAFT  CO.  Dcpl  501  <R  H.  Fayman),  Si  Louis  MO 
MEXICO  R  Cardenas 

MOHIl  PIPE  l  INI  CO  DAI  l  AS.  I  X  MGR  OE  ENGR  (NOACK) 

MCI  SI  R.  RUTI  EDGE,  WENTWORTH  AND  JOHNSI  ON  NEW  YORK  (RICH  ARDS) 

NEW  Z.EAI  AND  New  Zealand  Concrete  Research  Assoc  (Librarian).  Porirua 
NEWPORT  NEWS  SHIPBl  IK!  &  DRYIKICK  CO.  Newport  News  VA  (Tech.  Lib.) 

NORW  AY  DL  I  NORSKI.  VI  RI  TAS  (Korcn)  Oslo;  I.  l  oss.  Oslo;  J  Creed.  Ski;  Norwegian  Tech  llniv  (BrandUacg). 
Trondheim 

PACIFIC  MARINI-  TECHNOLOGY  Long  Beach.  CA  (Wagner) 

PORTLAND  CEMENT  ASSOC.  SKOKIK.  IHCORLF.Y;  Skokie  IL  (Rsch  &  Dev  Lab.  Lib.) 

PRESCON  CORP  TOW  SON.  MD  (KELLER) 

RANDCORP.  Santa  Monica  CA  (A.  I.aupa) 

RAYMOND  INTERNATIONAL  INC.  E  Colle  Soil  Tech  Dcpl.  Pennsaukcn.  NJ 
RIVERSIDE  CEMENT  CO  Riverside  CA  (W  Smith) 

SANDIA  LABORATORIES  Library  Div..  Livermore  CA 
SC  HO  PACK  ASSOC  SO  NORWALK.  CT  (SCHCPACK) 

SI  ATI  C  H  CORP.  MIAMI.  EL  (PERONI) 

SHELL  DEVEl  OPMENT  CO  Houston  TX  (E.  Doyle) 

she: 1 1  on  co.  Houston,  rx  (Marshall) 

SOUTH  AMERICA  N.  Nonet.  Valencia.  Venezuela 
SWEDEN  CicoTech  Inst;  VBB  (Library).  Stockholm 
TIDE  W  ATER  CONSTR.  CO  Norfolk  VA  (Eowlei ) 

TRW  SYSTEMS  Cl  EVELAND.  OH  (ENG.  I  IB  ):  REDONDO  BEACH.  CA  (DAI) 

UNI  TI  D  KINGDOM  Cement  &  Concrete  Assoc  W  exham  Springs,  Slough  Bucks;  Cement  &  Concrete  Assoc  (I  it. 
E\>.  Bucks;  D  New.  G.  Maunsell  Jk  Partners.  London:  I  ibrary,  Bristol;  R.  Browne.  Southall.  Middlesex;  Shaw  .V 
Hatton  (E  Hansen).  London;  I  aylor.  Woodrow  Conslr  (0I4P).  Southall.  Middlesex;  Unix,  of  Bristol  tR  Moigan). 
Bristol 

WAIT  BRIAN  ASSOC  INC.  Houston.  IX 

WES  TING  HOUSE  El  I  CERIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib.  Bryan);  Library.  Pittsburgh  PA 
WISS.  JANNEY.  ELSTNER.  &  ASStK'  Northbrook.  IL  (D.W.  Pfeifer) 

W(X)DWARl>Cl  YDE  CONSULTANTS  (A.  Harrigan)  San  Francisco;  PI  YMOUTH  MEETING  PA  (('ROSS.  Ill) 

Al  SMOOTS  Los  Angeles.  CA 
BARA.  JOHN  P  Lakewood.  CO 
BROWN.  ROBERT  University.  AL 
BULLOCK  La  Canada 
ERVIN.  DOUG  Belmont.  CA 
E  HEUZE  Boulder  CO 
LAYTON  Redmond.  WA 
CAPT  MURPHY  Sunnyvale.  CA 
R  E.  BESIER  Old  Saybrook  CT 
T.W.  MERMEl.  Washington  DC 
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